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The nature of the globule-coil transition of surface-confined polymers has been an issue of debate. 
Here this 2D coflapse transition is studied through a partially directed lattice model. In the general 
case of polymers with positive bending stiffness (A > 0), the collapse transition is first-order, it 
becomes second-order only in the limiting case of zero bending stiffness (A = 0). These analytical 
results are confirmed by Monte Carlo simulations. We also suggest some possible future experiments. 
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The collapse (globule-coil) transition of a self- 
attracting chain is one of the fundamental problems in 
polymer physics [H, 0, Q . Being deeply connected with 
biophysical problems such as DNA condensation, chro- 
matin organization, and protein folding, the collapse 
transition has also biological relevance. Generally speak- 
ing, this transition is caused by the competition be- 
tween monomer- monomer attraction and configurational 
entropy: Formation of contacts lowers the energy, but it 
requires monomers to be aligned and close to each other, 
thereby decreasing the polymer's degrees of configura- 
tional freedom. When the monomer- monomer attraction 
is much stronger than thermal energy, the polymer takes 
globular conformations to maximize the number of con- 
tacts. At the other limit of high temperatures, the poly- 
mer is in completely disordered coil states with maximal 
entropy. At zero external force, the collapse transition 
occurs at a so-called ^-temperature [l| . 

The exact nature of the collapse transition, however, 
is not yet completely settled, despite decades of exten- 
sive efforts [H, H, Hj]. This transition can be studied by 
mapping to the tricriticalpoint of the 4 -</> 6 0{n) field 
theory at the n~^0 limit [j, d, 0] , and is expected to be 
second-order in 2D and beyond (see, e.g., [7J]); further- 
more, exponents of the temperature-driven collapse and 
adsorption of a 2D polymer grafted on a linear bound- 
ary was obtained by analytical calculations [1, Q . How- 
ever, experimentally observed collapse transitions of both 
relatively flexible [l(| and semiflexible [ll[ 3D polymers 
are more like first-order transitions. Recent Monte Carlo 
(MC) simulation [13] is also in favor of the latter inter- 
pretation. 

With the advancement in single-molecule force manip- 
ulation methods, now the globule-coil transition can also 
be induced by external stretching 13j, [14J . According to 
mean-field theory the force-induced collapse transi- 
tion is first-order in all dimensions. This claim is con- 



firmed by MC simulation in 3D [16[, while in 2D the 
transition is argued to be second-order by MC simula- 
tion and scaling analysis 17. 18j|. 

In order to understand the collapse transition more 
deeply, here we investigate a 2D partially directed lat- 
tice model of a polymer chain which is exactly solvable 
(Fig. [I}. As the free energy density of this model system 
can be calculated precisely, we are able to draw definite 
conclusions concerning the nature of the collapse transi- 
tion. We find that the collapse transition is second-order 
for a polymer with exactly zero bending stiffness, while 
it changes to be first-order when the bending stiffness A 
is non-zero. Therefore bending energies have a dramatic 
effect on the cooperativity of the 2D globule-coil transi- 
tion. We also substantiate our predictions by performing 
extensive MC simulations. It becomes experimentally 
feasible to confine polymers to 2D mobile surfaces (see. 
e.g., (HI), therefore, future experiments will be able to 
verify the theoretical predictions. The present work is 
also directly applicable to the biology-related problem of 
protein /3-sheet unfolding. 

The model. — The 2D partially directed polymer of N 
identical units on a square lattice is shown in Fig. [T] The 
length of the bond connecting two consecutive monomers 
i and i+l is fixed to eto, while the direction of bond in — zo 
is prohibited. If any two monomers i and i+m (m>3) 
occupy nearest neighboring lattice sites, an attractive en- 
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FIG. 1: The 2D lattice polymer model with contacting inter- 
action, bending stiffness, and external stretching. The arrow 
shows the zo direction. 
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ergy of magnitude e is gained. Usually real polymers are 
semiflexible, we associate an energy penalty of magnitude 
A to each local direction change of the chain [2(| [2l| . For 
this model, only the special case of zero bending stiffness 
has been studied analytically 18, [12, Hj|. The results 



of the present paper suggest that, in the general case of 
A>0 the properties of the polymer will be dramatically 
different from these earlier predictions HI, 22> 23 1 . 

Partition functions. — To calculate the free energy den- 
sity of the polymer, a given configuration of the 2D chain 
is divided into a linear sequence of /3-sheet segments and 
coil segments [24 |. A /3-sheet segment is defined as a 
folded segment of np>2 consecutive columns, in which 
contacting interactions exist between any two adjacent 
ones. Two consecutive /3-sheet segments are separated 
by a coil segment, which is a segment of n c >0 consecu- 
tive columns in which all monomers are free of contacts. 
For example, the configuration shown in Fig. [T] has two 
/3-sheets and two coils. After having made such a dis- 
tinction between /3-sheets and coils, we proceed by first 
calculating the partition functions of /3-sheets and coils 
separately. 

Under the action of an external force /, the energy of 
a /3-sheet of np columns is 
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where lj is the number of monomers in the j-th column 
of the /3-sheet, and v(Z 3 -,Zj+i) = min(Zj-, Zj+i) — 1. The 
partition function of a /3-sheet with n>A monomers is 

[n/2] m-1 
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where p — e^ a °/ T , s = e A / T j a — e e / T , and T is the tem- 
perature. It is easier to calculate the generating function 
Gp(t^) of the partition function Zp(n) than to calculate 
Zp(n) directly. After some simple matrix operations 25 1 
we find that 



Gp(C) = J2((/a) n Zp(r 
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InEq. Q, Xj = (C/a) (j+1)/2 ; Ai(C) > A 2 (C) > ... are the 
eigenvalues of a L x L real-symmetric matrix A(C) with 
elements Ay(C) = C 1+(i+3 ' )/2 /o l<_il/2 = 1,2,..., L); 
and the orthogonal matrix contains the eigenvec- 

tors of matrix A(^). The parameter L should be infinity. 
When £ > 1, Ai(() = +oo and consequently Gp(() is 
not properly defined. When £ < 1, all the eigenvalues 
of matrix A(£) are finite, and the value of Gp(() can be 



calculated by Eq. (|4|). The point £ = 1 is therefore a 
singular point of Gp(Q, where it attains a finite value as 
long as Ai(l) < a/ (ps 2 ). At this point, the first derivative 
of Gp(Q with respect to £ is [25J 

L 

G'p(l) = dGp(C)/dC| C =i = lim Yl Mn) + 0(1) , (5) 

n=l 

where y(n) = (J2f.k=i A njA kj x k \j / (1 -ps 2 \j/a)) 2 . 



The configurational energy of a coil segment is 
E c = -n c fa + m c A 



(6) 



where n c and m c are, respectively, the total number of 
columns and the total number of bends in the configura- 
tion. In a coil segment the value of m c strongly depdends 
on the configurations. To calculate the partition function 
Z c (n) of a coil segment of n monomers, one needs to dis- 
tinguish among four different boundary conditions [25j |. 
The generating function G C (C) of Z c (n) in this case is 
expressed as 



G C (C) = YiQ/afZM 
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The divergence radius of G c (£) is easy to obtain from 
Eq. ©. G C (Q approaches +00 as C approaches this di- 
vergence radius. 

The collapse transition. — Since every configuration of 
the polymer is of the form . . .-/3-c-/3-c. . ., the generating 
function G(£) of the whole polymer's partition function 
Z(N) can readily be written down 3, 25j ]: 
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Z(iV) is related to the free energy density g(f, T) of the 
polymer by Z(N) = exp (-Nj3g(f,T)). Therefore, in the 
thermodynamic limit N — > 00, the free energy density is 



<?(/,T) = -e + TlnCo , 



(10) 



where Co is the smallest positive root in the range < 
Co < 1 of the following equation 



G /3 (Co)G c (Co) = 1 



(11) 



If Eq. (fTT|) has no root in the range of < Co < 1, then 
Co = I- The polymer's relative extension £ along zq is 
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A similar expression can also be obtained for the den- 
sity of contacts. For the simplicity we take the relative 
extension £ as our order parameter. 
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When both temperature T and external force / are 
sufficiently low, Eq. (fTTj) has no root in < Co < 1- Then 
g(f, T) = — e. The polymer is in the /3-sheet phase, and 
I = 0. There is no significant entropic contribution to 
the free energy density. As the temperature or the exter- 
nal force is elevated to certain point such that Eq. (fTl) is 
satisfied exactly at Co — 1, a phase transition occurs. At 
this point, the /9-sheet phase change to the extended coil 
phase. According to the work of Lifson [511] , this globule- 
coil transition is first-order if the value of G'p ( 1) in Eq. (J5J 
is finite at the phase transition point; otherwise, it is a 
continuous phase transition. At the collapse phase tran- 
sition point, high-precision numerical calculations reveal 
that y(n) of Eq. ((5]) decays exponentially with n as long 
as A > [5^ . In the vicinity of A=0 we have 



3 -0.88nA/e 



(A « e) 
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Therefore, at the phase transition point G'n(l) is finite as 
long as A > 0. In other words, the collapse phase tran- 
sition is first-order for polymers with bending stiffness 
A > 0; it is second-order for polymers with exactly zero 
bending stiffness. It is remarkable that the cooperativity 
of the collapse transition can be changed by adding just 
a small bending stiffness. 
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FIG. 2: (Color online) Response behavior of the polymer. 
(A) temperature-extension curve at force / = 0; (B) Force- 
extension curve at temperature T = 0.590928e. Solid lines are 
analytical predictions and symbols are MC simulation results 
for different system sizes. 



Results for A = 0. We perform analytical calcula- 
tions as well as Monte Carlo simulations on the model 
system. The £-T relation at / = is shown in Fig. [2^, 
and the £-f relation at T = 0.591e is shown in Fig. [2t>- 
When there is no external force, a second-order globule- 
coil phase transition occurs at T gc (0) = 0.8205e. The 
force-induced globule-coil transition at constant temper- 
ature is also second-order. Our findings are in accordance 
with Refs. 3, HI . In the vicinity of the collapse transi- 
tion, the relative extension of the extended coil phase can 
be expressed in a scaling form in terms of temperature 
change ST or force change Sf: 



£{T sc (f) + ST) 
£(f gc (T) + Sf) 



£(T gc (/))+ Cl (ST)-* , (14) 
£(f gc (T))+c 2 (6f)T , (15) 



where T gc (/) and f gc (T) are, respectively, the transition 
temperature at fixed force and the transition force at 
fixed temperature; 71 and 72 are two scaling exponents; 
and Ci,C2 are two numerical constants. In the case of 
A = 0, we have £(T gc {f)) = £(f gc (T)) = 0; and we 



find that the scaling exponent in Eq. pT|) is independent 
of force with 71 = 1/2 (see Fig. [3]), which is consistent 
with a second-order phase transition. The scaling rela- 
tion Eq. (|15j) with 72 = 1/2 is also confirmed by numer- 
ical calculation (data not shown). 

Results for A > 0. In this case, the collapse transi- 
tion becomes first-order. For example at A = 0.5e and 
/ = 0, the relative extension jumps from zero to £ = 0.193 
at the transition temperature T gc (0) = 1.209e (Fig. [2^). 
Such a large jump is also observed in the force-induced 
transition (Fig. [2Jd). A non-zero bending stiffness there- 
fore is able to dramatically enhance the cooperativity of 
the globule-coil phase transition. This may be partially 
understood in the following way. A positive bending en- 
ergy significantly decreases the configurational entropy 
of a coil segment. Consequently the globule-coil transi- 
tion will occur at higher temperature and higher force, 
and once the polymer is unfolded it favors those highly 
elongated configurations which have fewer bends. 

As is expected for a first-order phase transition, in the 
case of A > 0, the scaling exponents in the scaling forms 
and (H5J are 71 = 72 = 1 (Fig. g). 

Conclusion. — In this work we studied the collapse tran- 
sition of a 2D partially directed lattice polymer model 
using the Lifson approach [5i|. We devided each con- 
figuration of the polymer into /3-sheet segments and coil 
segments and then calculated their partition functions 
separately. This approach enables us to precisely cal- 
culate the free energy density of the system. Previous 
tricritical field theory [55|, [53] and numerical scaling anal- 
ysis [HI] predicted that the collapse transition of a flex- 
ible (A=0) partially directed lattice polymer is second- 
order. We have confirmed this prediction. Our calcula- 
tions further revealed that, the susceptibility \ — d£/df 

— 1/2 

of the polymer diverges as x oc (/ — f gc (T)) or 
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loglO(8T/e) 



a mobile lipid bilayer [19j |. Polymer configurations can 
be recorded in real-time and manipulated by controlling 
temperature or external force. 
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FIG. 3: The scaling relation between extension increment Si 
[= l(T gc (f) + ST) - e{T gc (f))] and ST. The upper and lower 
dotted line has slope 1/2 and 1, respectively. The dimension 
of the matrix A(£) is set to L = 5000. 
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(T — T gc (f)) as / or T approaches the phase 
transition value from above. The most important result 
of the present paper is that, a non-zero bending stiff- 
ness of the polymer changes the nature of the collapse 
transition from being second-order to being first-order. 
This conclusion is in agreement with an earlier exact enu- 
meration study [2^1 . It is well known that, in the 2D 
Ising model the paramagnet-ferromagnet phase transi- 
tion changes from being second-order to being first-order 
under the action of a non-zero external magnetic field. 
It it is interesting that in the 2D polymer system such a 
qualitative change is caused not by an external field (such 
as the external force), but by an internal (microscopic) 
parameter, the bending stiffness A. 

At this point, it remains to be seen whether the above 
conclusion also holds if we remove the constraint of par- 
tial directedness. In the case of DNA thermal denatura- 
tion, including long-range excluded-volume interactions 
between different bubbles changes the order of the phase 
transition 26]. When additional long-range excluded- 
volume interactions are included, it may not be surprising 
if the temperature-induced 2D collapse transition will be 
first-order only when the bending energy is much larger 
than the contacting energy e 2l|. On the other hand, 
the force-induced 2D collapse transition may have differ- 
ent behavior. Mean-field work of Ref. [3] suggested that 
in the globule-coil transition, when a force is approached 
from below, a collapsed 2D polymer takes shape of an 
elongated elipse (i.e., partially directed). The mathe- 
matical constraint of partial directedness of the present 
paper may not qualitatively change the cooperativity of 
the force-induced 2D collapse transition. Additional MC 
simulation work as well as exact numerical enumeration 
study are needed to clarify this point. (When A = 0.3e, 
our prelimilary MC results hinted at a first-order phase 
transition, with f c 0.9e at T = 0.6e.) It is also highly 
desirable to perform real 2D polymer collapse experi- 
ments. For example, a flexible or semiflexible biopoly- 
mer [such as DNA and poly (Glycine)] can be attached to 
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